On estimating the basic reproduction number in distinct stages of a contagious disease spreading

In epidemiology, the basic reproduction number R-0 is usually defined as the average number of new infections caused by a single infective individual introduced into a completely susceptible population. According to this definition. R-0 is related to the initial stage of the spreading of a contagious disease. However, from epidemiological models based on ordinary differential equations (ODE), R-0 is commonly derived from a linear stability analysis and interpreted as a bifurcation parameter: typically, when R-0 >1, the contagious disease tends to persist in the population because the endemic stationary solution is asymptotically stable: when R-0 <1, the corresponding pathogen tends to naturally disappear because the disease-free stationary solution is asymptotically stable. Here we intend to answer the following question: Do these two different approaches for calculating R-0 give the same numerical values? In other words, is the number of secondary infections caused by a unique sick individual equal to the threshold obtained from stability analysis of steady states of ODE? For finding the answer, we use a susceptibleinfective-recovered (SIR) model described in terms of ODE and also in terms of a probabilistic cellular automaton (PCA), where each individual (corresponding to a cell of the PCA lattice) is connected to others by a random network favoring local contacts. The values of R-0 obtained from both approaches are compared, showing good agreement. (C) 2012 Elsevier B.V. All rights reserved.

Vectorial capacity, basic reproduction number, force of infection and all that: formal notation to complete and adjust their classical concepts and equations

A dimensional analysis of the classical equations related to the dynamics of vector-borne infections is presented. It is provided a formal notation to complete the expressions for the Ross' threshold theorem, the Macdonald's basic reproduction "rate" and sporozoite "rate", Garret-Jones' vectorial capacity and Dietz-Molineaux-Thomas' force of infection. The analysis was intended to provide a formal notation that complete the classical equations proposed by these authors.

Estimating the Reproduction Number of Ebola Virus (EBOV) During the 2014 Outbreak in West Africa

The 2014 Ebola virus (EBOV) outbreak in West Africa is the largest outbreak of the genus Ebolavirus to date. To better understand the spread of infection in the affected countries, it is crucial to know the number of secondary cases generated by an infected index case in the absence and presence of control measures, i.e., the basic and effective reproduction number. In this study, I describe the EBOV epidemic using an SEIR (susceptible-exposed-infectious-recovered) model and fit the model to the most recent reported data of infected cases and deaths in Guinea, Sierra Leone and Liberia. The maximum likelihood estimates of the basic reproduction number are 1.51 (95% confidence interval [CI]: 1.50-1.52) for Guinea, 2.53 (95% CI: 2.41-2.67) for Sierra Leone and 1.59 (95% CI: 1.57-1.60) for Liberia. The model indicates that in Guinea and Sierra Leone the effective reproduction number might have dropped to around unity by the end of May and July 2014, respectively. In Liberia, however, the model estimates no decline in the effective reproduction number by end-August 2014. This suggests that control efforts in Liberia need to be improved substantially in order to stop the current outbreak.

Rapid drop in the reproduction number during the Ebola outbreak in the Democratic Republic of Congo.

The Democratic Republic of Congo (DRC) experienced a confined rural outbreak of Ebola virus disease (EVD) with 69 reported cases from July to October 2014. Understanding the transmission dynamics during the outbreak can provide important information for anticipating and controlling future EVD epidemics. I fitted an EVD transmission model to previously published data of this outbreak and estimated the basic reproduction number R 0 = 5.2 (95% CI [4.0-6.7]). The model suggests that the net reproduction number Rt fell below unity 28 days (95% CI [25-34] days) after the onset of symptoms in the index case. This study adds to previous epidemiological descriptions of the 2014 EVD outbreak in DRC, and is consistent with the notion that a rapid implementation of control interventions helped reduce further spread.

A dimensionless number for understanding the evolutionary dynamics of antigenically variable RNA viruses.

Antigenically variable RNA viruses are significant contributors to the burden of infectious disease worldwide. One reason for their ubiquity is their ability to escape herd immunity through rapid antigenic evolution and thereby to reinfect previously infected hosts. However, the ways in which these viruses evolve antigenically are highly diverse. Some have only limited diversity in the long-run, with every emergence of a new antigenic variant coupled with a replacement of the older variant. Other viruses rapidly accumulate antigenic diversity over time. Others still exhibit dynamics that can be considered evolutionary intermediates between these two extremes. Here, we present a theoretical framework that aims to understand these differences in evolutionary patterns by considering a virus's epidemiological dynamics in a given host population. Our framework, based on a dimensionless number, probabilistically anticipates patterns of viral antigenic diversification and thereby quantifies a virus's evolutionary potential. It is therefore similar in spirit to the basic reproduction number, the well-known dimensionless number which quantifies a pathogen's reproductive potential. We further outline how our theoretical framework can be applied to empirical viral systems, using influenza A/H3N2 as a case study. We end with predictions of our framework and work that remains to be done to further integrate viral evolutionary dynamics with disease ecology.

On the time to reach a critical number of infections in epidemic models with infective and susceptible immigrants.

In this paper we examine the time T to reach a critical number K0 of infections during an outbreak in an epidemic model with infective and susceptible immigrants. The underlying process X, which was first introduced by Ridler-Rowe (1967), is related to recurrent diseases and it appears to be analytically intractable. We present an approximating model inspired from the use of extreme values, and we derive formulae for the Laplace-Stieltjes transform of T and its moments, which are evaluated by using an iterative procedure. Numerical examples are presented to illustrate the effects of the contact and removal rates on the expected values of T and the threshold K0, when the initial time instant corresponds to an invasion time. We also study the exact reproduction number Rexact,0 and the population transmission number Rp, which are random versions of the basic reproduction number R0.

A simple model for the establishment of tick-borne pathogens of Ixodes scapularis : a global sensitivity analysis of R0

The basic reproduction number of a pathogen, R 0, determines whether a pathogen will spread (R0>1R 0>1), when introduced into a fully susceptible population or fade out (R0<1R 0<1), because infected hosts do not, on average, replace themselves. In this paper we develop a simple mechanistic model for the basic reproduction number for a group of tick-borne pathogens that wholly, or almost wholly, depend on horizontal transmission to and from vertebrate hosts. This group includes the causative agent of Lyme disease, Borrelia burgdorferi, and the causative agent of human babesiosis, Babesia microti, for which transmission between co-feeding ticks and vertical transmission from adult female ticks are both negligible. The model has only 19 parameters, all of which have a clear biological interpretation and can be estimated from laboratory or field data. The model takes into account the transmission efficiency from the vertebrate host as a function of the days since infection, in part because of the potential for this dynamic to interact with tick phenology, which is also included in the model. This sets the model apart from previous, similar models for R0 for tick-borne pathogens. We then define parameter ranges for the 19 parameters using estimates from the literature, as well as laboratory and field data, and perform a global sensitivity analysis of the model. This enables us to rank the importance of the parameters in terms of their contribution to the observed variation in R0. We conclude that the transmission efficiency from the vertebrate host to Ixodes scapularis ticks, the survival rate of Ixodes scapularis from fed larva to feeding nymph, and the fraction of nymphs finding a competent host, are the most influential factors for R0. This contrasts with other vector borne pathogens where it is usually the abundance of the vector or host, or the vector-to-host ratio, that determine conditions for emergence. These results are a step towards a better understanding of the geographical expansion of currently emerging horizontally transmitted tick-borne pathogens such as Babesia microti, as well as providing a firmer scientific basis for targeted use of acaricide or the application of wildlife vaccines that are currently in development.

A Generalist Predator Regulating Spread of a Wildlife Disease: Exploring Two Infection Transmission Scenarios

Ecoepidemiology is a well-developed branch of theoretical ecology, which explores interplay between the trophic interactions and the disease spread. In most ecoepidemiological models, however, the authors assume the predator to be a specialist, which consumes only a single prey species. In few existing papers, in which the predator was suggested to be a generalist, the alternative food supply was always considered to be constant. This is obviously a simplification of reality, since predators can often choose between a number of different prey. Consumption of these alternative prey can dramatically change their densities and strongly influence the model predictions. In this paper, we try to bridge the gap and explore a generic ecoepidemiological system with a generalist predator, where the densities of all prey are dynamical variables. The model consists of two prey species, one of which is subject to an infectious disease, and a predator, which consumes both prey species. We investigate two main scenarios of infection transmission mode: (i) the disease transmission rate is predator independent and (ii) the transmission rate is a function of predator density. For both scenarios we fulfil an extensive bifurcation analysis. We show that including a second dynamical prey in the system can drastically change the dynamics of the single prey case. In particular, the presence of a second prey impedes disease spread by decreasing the basic reproduction number and can result in a substantial drop of the disease prevalence. We demonstrate that with efficient consumption of the second prey species by the predator, the predator-dependent disease transmission can not destabilize interactions, as in the case with a specialist predator. Interestingly, even if the population of the second prey eventually vanishes and only one prey species finally remains, the system with two prey species may exhibit different properties to those of the single prey system.

Control dynamics of severe acute respiratory syndrome transmission

Severe acute respiratory syndrome (SARS) is a serious disease with many puzzling features. We present a simple, dynamic model to assess the epidemic potential of SARS and the effectiveness of control measures. With this model, we analysed the SARS epidemic data in Beijing. The data fitting gives the basic case reproduction number of 2.16 leading to the outbreak, and the variation of the effective reproduction number reflecting the control effect. Noticeably, our study shows that the response time and the strength of control measures have significant effects on the scale of the outbreak and the lasting time of the epidemic.

Investigating the persistence of tick-borne pathogens via the R-0 model

In the epidemiology of infectious diseases, the basic reproduction number, R-0, has a number of important applications, most notably it can be used to predict whether a pathogen is likely to become established, or persist, in a given area. We used the R-0 model to investigate the persistence of 3 tick-borne pathogens; Babesia microti, Anaplasma phagocytophilum and Borrelia burgdorferi sensu lato in an Apodemus sylvaticus-Ixodes ricinus system. The persistence of these pathogens was also determined empirically by screening questing ticks and wood mice by PCR. All 3 pathogens behaved differently in response to changes in the proportion of transmission hosts on which I. ricinus fed, the efficiency of transmission between the host and ticks and the abundance of larval and nymphal ticks found on small mammals. Empirical data supported theoretical predictions of the R-0 model. The transmission pathway employed and the duration of systemic infection were also identified as important factors responsible for establishment or persistence of tick-borne pathogens in a given tick-host system. The current study demonstrates how the R-0 model can be put to practical use to investigate factors affecting tick-borne pathogen persistence, which has important implications for animal and human health worldwide.

Optimal control and numerical optimization applied to epidemiological models

A relação entre a epidemiologia, a modelação matemática e as ferramentas computacionais permite construir e testar teorias sobre o desenvolvimento e combate de uma doença. Esta tese tem como motivação o estudo de modelos epidemiológicos aplicados a doenças infeciosas numa perspetiva de Controlo Ótimo, dando particular relevância ao Dengue. Sendo uma doença tropical e subtropical transmitida por mosquitos, afecta cerca de 100 milhões de pessoas por ano, e é considerada pela Organização Mundial de Saúde como uma grande preocupação para a saúde pública. Os modelos matemáticos desenvolvidos e testados neste trabalho, baseiam-se em equações diferenciais ordinárias que descrevem a dinâmica subjacente à doença nomeadamente a interação entre humanos e mosquitos. É feito um estudo analítico dos mesmos relativamente aos pontos de equilíbrio, sua estabilidade e número básico de reprodução. A propagação do Dengue pode ser atenuada através de medidas de controlo do vetor transmissor, tais como o uso de inseticidas específicos e campanhas educacionais. Como o desenvolvimento de uma potencial vacina tem sido uma aposta mundial recente, são propostos modelos baseados na simulação de um hipotético processo de vacinação numa população. Tendo por base a teoria de Controlo Ótimo, são analisadas as estratégias ótimas para o uso destes controlos e respetivas repercussões na redução/erradicação da doença aquando de um surto na população, considerando uma abordagem bioeconómica. Os problemas formulados são resolvidos numericamente usando métodos diretos e indiretos. Os primeiros discretizam o problema reformulando-o num problema de optimização não linear. Os métodos indiretos usam o Princípio do Máximo de Pontryagin como condição necessária para encontrar a curva ótima para o respetivo controlo. Nestas duas estratégias utilizam-se vários pacotes de software numérico. Ao longo deste trabalho, houve sempre um compromisso entre o realismo dos modelos epidemiológicos e a sua tratabilidade em termos matemáticos.

A mathematical model for the dynamics of Banana Xanthomonas Wilt with vertical transmission and inflorescence infection

A mathematical model for Banana Xanthomonas Wilt (BXW) spread by insect is presented. The model incorporates inflorescence infection and vertical transmission from the mother corm to attached suckers, but not tool-based transmission by humans. Expressions for the basic reproduction number R0 are obtained and it is verified that disease persists, at a unique endemic level, when R0 > 1. From sensitivity analysis, inflorescence infection rate and roguing rate were the parameters with most influence on disease persistence and equilibrium level. Vertical transmission parameters had less effect on persistence threshold values. Parameters were approximately estimated from field data. The model indicates that single stem removal is a feasible approach to eradication if spread is mainly via inflorescence infection. This requires continuous surveillance and debudding such that a 50% reduction in inflorescence infection and 2–3 weeks interval of surveillance would eventually lead to full recovery of banana plantations and hence improved production.

A mathematical model for the spread of streptococcus pneumoniae with transmission due to sequence type

This paper discusses a simple mathematical model to describe the spread of Streptococcus pneumoniae. We suppose that the transmission of the bacterium is determined by multi-locus sequence type. The model includes vaccination and is designed to examine what happens in a vaccinated population if MLSTs can exist as both vaccine and non vaccine serotypes with capsular switching possible from the former to the latter. We start off with a discussion of Streptococcus pneumoniae and a review of previous work. We propose a simple mathematical model with two sequence types and then perform an equilibrium and (global) stability analysis on the model. We show that in general there are only three equilibria, the carriage-free equilibrium and two carriage equilibria. If the effective reproduction number Re is less than or equal to one, then the carriage will die out. If Re > 1, then the carriage will tend to the carriage equilibrium corresponding to the multi-locus sequence type with the largest transmission parameter. In the case where both multi-locus sequence types have the same transmission parameter then there is a line of carriage equilibria. Provided that carriage is initially present then as time progresses the carriage will approach a point on this line. The results generalize to many competing sequence types. Simulations with realistic parameter values confirm the analytical results.

Biodiversity Can Help Prevent Malaria Outbreaks in Tropical Forests

Background: Plasmodium vivax is a widely distributed, neglected parasite that can cause malaria and death in tropical areas. It is associated with an estimated 80-300 million cases of malaria worldwide. Brazilian tropical rain forests encompass host- and vector-rich communities, in which two hypothetical mechanisms could play a role in the dynamics of malaria transmission. The first mechanism is the dilution effect caused by presence of wild warm-blooded animals, which can act as dead-end hosts to Plasmodium parasites. The second is diffuse mosquito vector competition, in which vector and non-vector mosquito species compete for blood feeding upon a defensive host. Considering that the World Health Organization Malaria Eradication Research Agenda calls for novel strategies to eliminate malaria transmission locally, we used mathematical modeling to assess those two mechanisms in a pristine tropical rain forest, where the primary vector is present but malaria is absent. Methodology/Principal Findings: The Ross-Macdonald model and a biodiversity-oriented model were parameterized using newly collected data and data from the literature. The basic reproduction number (R0) estimated employing Ross-Macdonald model indicated that malaria cases occur in the study location. However, no malaria cases have been reported since 1980. In contrast, the biodiversity-oriented model corroborated the absence of malaria transmission. In addition, the diffuse competition mechanism was negatively correlated with the risk of malaria transmission, which suggests a protective effect provided by the forest ecosystem. There is a non-linear, unimodal correlation between the mechanism of dead-end transmission of parasites and the risk of malaria transmission, suggesting a protective effect only under certain circumstances (e.g., a high abundance of wild warm-blooded animals). Conclusions/Significance: To achieve biological conservation and to eliminate Plasmodium parasites in human populations, the World Health Organization Malaria Eradication Research Agenda should take biodiversity issues into consideration. © 2013 Laporta et al.

Uma abordagem probabilística do número de reprodução básica em modelos epidemiológicos com aplicação na ferrugem do eucalipto

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)